Curvature dependence of the electrolytic liquid-liquid interfacial tension 
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The interfacial tension of a liquid droplet surrounded by another hquid in the presence of micro- 
scopic ions is studied as a function of the droplet radius. An analytical expression for the interfacial 
tension is obtained within a linear Poisson-Boltzmann theory and compared with numerical re- 
sults from non-linear Poisson-Boltzmann theory. The excess liquid-liquid interfacial tension with 
respect to the pure, salt-free liquid-liquid interfacial tension is found to decompose into a curvature- 
independent part due to short-ranged interfacial effects and a curvature-dependent electrostatic 
contribution. Several curvature-dependent regimes of different scaling of the electrostatic excess in- 
terfacial tension are identified. Symmetry relations of the interfacial tension upon swapping droplet 
and bulk liquid are found to hold in the low-curvature limit, which, e.g., lead to a sign change of the 
excess Tolman length. For some systems a low-curvature expansion up to second order turns out 
to be applicable if and only if the droplet size exceeds the Debye screening length in the droplet, 
independent of the Debye length in the bulk. 

PACS numbers; 68.05.-n, 68.03. Cd, 82.45.Gj 



I. INTRODUCTION 

Common wisdom in emulsion science tells that, in 
order to kinetically stabilize an emulsion of water and 
oil, say, surfactants are needed in order to decrease 
the interfacial tension thereby decreasing the thermody- 
namic force causing droplet coalescence This pic- 

\ ture has been upset by Leunissen et al. who showed ex- 
perimentally that in certain additive-free water-oil mix- 
tures micron-sized water droplets in oil may be stabi- 
lized electrostatically by absorbing ions present in the 
system 0, Q . Several aspects of these experiments such 
as the proposed chargingof the water droplets due to 
an unequal partitioning [2, [^, Q and the formation of 
a colloidal crystal of water droplets [3l can be under- 
stood theoretically within a simple Poisson-Boltzmann 
model. However, the rather unimodal size distribution of 
the water droplets in the above-mentioned experiments 
has not been explained so far. A similar observation has 
been made by Sacanna et al. who found experimental in- 
dications of the existence of thermodynamically favored 

' droplet radii in certain emulsions stabilized by nano-sized 
colloids [1| . A thermodynamically favored droplet radius 
requires a radius dependent water-oil interfacial tension 
because otherwise the global minimum of the free energy 
would be attained for one single macroscopic drop. One is 
thereby led to the problem of analyzing the liquid-liquid 
interfacial tension as a function of the droplet radius. 

The study of the curvature dependence of liquid- vapor 
surface tensions has been pioneered by Gibbs Q , Tolman 
[To| . and Kirkwood and Buff Tolman introduced 

a low-curvature expansion of the form 7(a)/7(oo) ~ 
l/(l + 2(5/a) ~ 1 — 2(5/a where a denotes the radius of 
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curvature, 7(a) is the surface tension of the curved sur- 
face, and 7(00) is its planar value. The parameter 5, 
which has the dimension of length, is called the Tolman 
length and it can be identified with the spatial distance 
between the Gibbs dividing surface and the surface of 
tension. In the last decades the concept of a curvature 
dependent liquid- vapor surface tension has been taken up 
within various studies on critical phenomena [l^ . inter- 
face elasticity [l^ , and nucleation [3, [l^ . 

However, whereas in all these investigations the droplet 
and the surrounding bulk were composed of the same 
substance, albeit in different phases, here a mixture of 
two different liquids and ions is studied. Moreover, only 
the excess interfacial tension due to the electrolyte is of 
interest here while the two liquids forming droplet and 
bulk merely act as external fields onto the ions. 

The present investigation is carried out within the 
spherical version of the model studied in Ref. [Si] f Sec. HI)) . 
As in Ref. H] linearization of the Poisson-Boltzmann 
equation offers the possibility of closed analytical expres- 
sions for the interfacial tension (Sec. Illip . In Sec. IIVI 
the approximative analytical expressions for the inter- 
facial tension will be shown to at least qualitatively, in 
many realistic cases even quantitatively, agree with the 
numerical results obtained within the full, non-linear the- 
ory. The main conclusion will be that it is precisely the 
electrostatic contribution to the interfacial tension that 
brings about a curvature dependence, which, however, is 
usually insignificant to serve as an explanation for uni- 
modal radius distributions in the emulsions by Leunissen 
et al. mentioned above (Sec.|V|. On the other hand, the 
curvature-dependent electrostatic contribution to the in- 
terfacial tension can be expected to increase considerably 
in magnitude if highly charged colloids instead of mono- 
valent ions are present. Under these conditions, however, 
the approximations made in the present work are not a 
priori justified, and it is left for future studies to inves- 
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tigate the influence of valency on the qualitative picture 
to be drawn here, which corresponds to the low-valency 
limit. 

II. MODEL AND FORMALISM 

In the following dimensionless quantities are expessed 
in units of the thermal energy k^T, the elementary charge 

e, and the vacuum Bjerrum length £ — with 

the permeability of the vacuum £vac • Dimensionful quan- 
tities are denoted by the same symbol as the correspond- 
ing dimensionless quantities. 

Consider a liquid spherical droplet of radius a and rel- 
ative dielectric constant Sd surrounded by bulk liquid of 
relative dielectric constant Eh- Due to the spherical sym- 
metry of the setting the only relevant positional vari- 
able is the distance r S [0, oo) from the droplet center. 
Monovalent cations (-1- ions) and anions (— ions) are dis- 
tributed in both liquids. The difference in solvation free 
energy of a ± ion in the droplet with respect to the bulk 
liquid is denoted by f± , which, within the Born approx- 
imation M , can be estimated by f± = ( ) 

2a±\ed e\,l 

with the ion radius a±. As in Ref. 5] all interfacial ef- 
fects due to, e.g., smooth interfaces, finite ion size, van 
der Waals forces, and image charges, which are short 
ranged as compared to the electrostatic potential, are ac- 
counted for by introducing solvent-induced ion potentials 
^±(^) — /±0(a + s — 7") with 8 the Heaviside function. 
Note that the parameter s, which describes the radial 
offset of the solvent induced ion potentials V± with re- 
spect to the dielectric interface at r = a, can be positive 
and negative, depending on the net effect of the above- 
mentioned interfacial effects. More detailed representa- 
tions of the interfacial effects are possible at the expense 
of more phenomenological parameters [13, EM [3, but 
for the sake of convenience and because handy analytical 
expressions are desired the present most simple choice is 
made here. 

A convenient approach to calculate the interfacial ten- 
sion of the system under consideration is to first deter- 
mine the equilibrium ion number densi ty p rofiles q± by 
means of density functional theory [lO, l2ll . [2^ and then 
to infer the interfacial tension from inserting these equi- 
librium profiles into the grand potential density func- 
tional. Poisson-Boltzmann theory corresponds to the 
mean-field grand potential density functional 

"=±0 

+ Vo.{r) + ^cP{r,[g±])) (1) 

with fj,a the chemical potential of a ions and (f>{r, [g±]) 
the electrostatic potential functional at radius r, which 



fulfills the Poisson equation 

l(£(r)r2,/,'(r, [g±]))' = -in ^ ag^ir) (2) 

Q = ± 

subject to the boundary conditions (f>'{r — 0) = and 
(j)(r = oo) = 0, where a prime denotes a derivative with 
respect to r and e(r) :— ed6(a — r) + ehQ{r — a). The 
electrostatic potential is a continuous function of r, and 
at the dielectric interface (r = a) the radial component of 
the dielectric displacement is continuous: ed4''{i' ^ — 
£h<f)'{r \ a). 

Minimizing the density functional in Eq. ([1]) gives rise 
to the Euler-Lagrange equations 

Qa{r) = exp(^a - y^ir) - a(f){r, [g±])). (3) 

Due to the local charge neutrality in the bulk liquid far 
away from the droplet {g+{r — oo) = g-{r = oo)) one in- 
fers fi-^- ~ /i_ =: fi. Upon introducing the reference den- 
sities gf exp(Ai) and gf ■.=^gf exp{-iU + f^)/2), 
the sharp- kink reference density profile g'^'^^{r,x) := 
g\^^Q{r — x) + g'^^Qix — r) with the discontinuity lo- 
cated at radius x, and the shifted electrostatic potential 
ijj{r) := (j){r) — (j)]jQ(a + s — r) with the Donnan potential 
:— (/_ — /+)/2, the Euler-Lagrange equation Q can 
be rewritten as 

Qair) = g^^H^^ a + s) exp{-aip{r)). (4) 

Inserting Eq. ^ into the Poisson equation ^ leads to 
the Poisson-Boltzmann equation 

'0"(r) = K(r)^ sinh('0(r)), r^a,a + s (5) 

with K{r) y^8TTg^'^^{r, a + s)/e{r) the Debye screen- 
ing factor. Given a solution tp, the interfacial tension 
with respect to the dielectric interface at r = a in excess 
to the pure, salt-free liquid-liquid interfacial tension be- 
tween the droplet and the bulk liquid is, after inserting 
Eq. (g]) into Eq. ([T|), determined by 

,,cx_ ^[g±]-"[g''°^(-.«)] 

^ 47ra2 

OO 

= E Jdrr^(^ga{r)-g'^'^'ir,a) 

+ ^go.ir)HrA8±]))- (6) 

As solutions of the non-linear Poisson-Boltzmann equa- 
tion ^ in the spherical geometry can be obtained only 
numerically, the same holds for the excess interfacial ten- 
sion 7°^^ in Eq. ^ . However, upon linearizing the Euler- 
Lagrange equation ([4]) and the Poisson-Boltzmann equa- 
tion ^ one obtains analytical expressions for the excess 
interfacial tension 7'''', which will be derived in the next 
section. 
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III. LINEARIZED THEORY 

For a sufficiently small Donnan potential \4>d\ < 
1 the Euler-Lagrange equations dH) and the Poisson- 
Boltzmann equation ^ can be linearized leading to 



and 



a 6\a 



with 



a{r) 



47r 



(7) 



'4)"{r) ^ K{ryTjj{r), r^a,a + s, (8) 



respectively. Inserting both expressions into Eq. ([6]) one 
obtains 



(9) 



(10) 



the charge enclosed by a sphere of radius r around the 
origin per sphere surface area. 

The linear Poisson-Boltzmann equation ([5]) is analyti- 
cally soluble which gives rise to an expression of the form 



a{a + s) 



s 

1 + - 
a 



gbgb 
2tt 



rcf 



-F{s/a,Ki,a,n,p), (11) 



where 



Kb := y^STTg^^V^b, n := v^Ed/eb, and p 



\J O'd^ / Qh^ ■ The full scaling function which is 
recorded in the appendix, appears somewhat lengthy but 
is straightforward to obtain in principle. 

However, since the effective interfacial width parame- 
ter s is usually very much smaller than the droplet ra- 
dius and the local Debye lengths, \s\ <C a,K(r)~^, the 
first argument of the scaling function F can, within an 
excellent approximation, be set to zero. Inserting x ~ Q 
into Eqs. (|A.1|) and (|A.2|) leads to 



np + 



n{p — n) /n 



F{0,y,n,p) = 



cxp 



2t/^) (np 



n{p + n) ('n\'^ 

y W 



i — / P\ f 1 — n"\ 

+ np-\ h exp - 2y— — 1 + np 

y \ n/ \ y ' 



(12) 



At this level of approximation Eqs. 
to 



= 2.,-f(l-p^)-^a(a) 



and (fTTjl reduce 
(13) 



and 



<T{a) - 4>D\j^-^F{0, K^a, n,p), (14) 

respectively. According to Eq. ([T^ . the droplet charge 
per droplet surface area cr(a) is (almost) independent of 
the interfacial width s. On the other hand, the excess 
interfacial tension 7°'' in Eq. comprises a contribu- 
tion describing the ion exclusion due to the short-ranged 
interfacial effects 



7i( 



■.= 2sgfil-p^), 



(15) 



which is (essentially) linear in s and (almost) indepen- 
dent of the droplet radius a, as well as an electrostatic 
contribution 



7o: 



(PD , ^ 



(16) 



which is (almost) independent of the effective interfacial 
width s. 



In order to understand the involved dependence of the 
scaling function i^(0, y, n, p) on y it is useful to investigate 
the asymptotic behavior for large and small values of y. 
li y ^ yi := n/p, the terms in Eq. (fT2|) proportional to 
the exponentials may be neglected such that 



np 



FiO,y:^y^,n,p):^—^G{y,n,p) (17) 
l + np 



with 



Giy,n,p) 



n \ / 1 
1-- 1 + - 

py^ ^ y 



1 



1 - 
(1 -I- np)y 



1 - n 



(18) 



Defining y^ := and y„ :— 1 and noting that 

l + np 



1 -n" 
\ + np 



e [—2/1^,2/3 ] one infers from Eq. pT)) the leading 



4 



O 



SI, 

sT 

ST 




2/2 

y [log scale] 



FIG. 1: Scaling function F{0,y,n,p) as a function of y for 
the relation <S J/2 ^ crossover positions (see 

main text) in a log-log plot. The asymptotic regimes I-IV 
corresponding to Eqs. (|19|l and (I20p are apparent. For ^ 
y2 regime III is absent, and for y^ S> y^ also regime II. 



order asymptotic behavior 



F(0,y » y^,n,p) < 



np 


1 


+ np 




np 


1 


+ 




np 


1 


— n^ 



(I) y » v'^ 

y^' (11) ya"" « 2/ < Vs 

(ni)y«y2^ 



(19) 

The three cases considered in Eq. ([TO)) are exhaustive 
and mutually exclusive for y 3> because 2/2 ii- 
max(yf , ). If y < y^^ , Eq. ((HI) leads to 



F(0,y <yf ,n,p) ~ —y. 



(20) 



Figure [T] displays F(0, y, n,p) for the case Vi •^y2 
, where all four asymptotic regimes I-IV of Eqs. 
and ((20)) are apparent. If y^ Vi ^ Us , however, 
regime III in Fig. [T] is absent, and a crossover between 
regimes II and IV takes place at y = yf . Moreover, 
if yf ^ yg regime II is also absent, and F{0,y,n,p) 
exhibits a single crossover at y = y^^ between regimes I 
and IV. 

It will turn out in the next section that the analyti- 
cal expressions based on the linearized theory derived in 
the present section agree qualitatively, in typical cases 
even quantitatively, with numerically calculated interfa- 
cial tensions within the non-linear theory. 



IV. DISCUSSION 

Here the closed analytical expressions obtained within 
the linearized Poisson-Boltzmann theory of the previ- 
ous section are discussed and compared with numerical 
results obtained within the non-linear theory based on 



-70 



-140 



■ 1 1 ' 

En = 5 








o/w^ 











FIG. 2: Electrostatic contribution to the excess interfacial 
tension 7^^ in mixtures of oil (eo G {5, 7.5, 10}) and water 
(ew = 80) as a function of the radius a of an oil droplet in 
water (O/W, ascending curves) and a water droplet in oil 
(W/0, descending curves) with ion radii a+ = 0.36 nm and 
a- = 0.30 nm as well as an ionic strength in water = 1 mM. 
The thin solid curves are calculated by means of the analytical 
expressions within the linear theory of Sec. IIIII whereas the 
thick dotted curves are obtained by numerically solving the 
non-linear Poisson-Boltzmann equation ([5]). Upon swapping 
oil and water (O/W <-> W/O) the slope of the curves at = 
(planar system) changes its sign. 
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(b) W/O 
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7.5 
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0.0545 



TABLE I: Quantities n, p, and (f)D as well as the crossover 
values y^ and j/^ (see Sees. [IT] and IIIip within the Born ap- 
proximaion for (a) O/W and (b) W/O systems with the oil 
dielectric constant £0 € {5, 7.5, 10} and ion radii 0+ = 0.36 nm 
and a_ = 0.30 nm. 



Eqs. (|l])-(in]). Some of the numerical data presented here 
have already been considered in Ref . . Throughout this 
section one of the liquids is water with dielectric constant 
£w = 80. Moreover, the largely arbitrary but representa- 
tive choice of ion radii a+ = 0.36 nm and a_ = 0.30 nm 
is made throughout. Given the dielectric constant of the 
second liquid, called "oil" , the parameters n, p, and <f)£) 
are known within the Born approximation (see Sec. |TT|. 
The cases of an oil droplet in water (O /W) and of a water 
droplet in oil (W/O) will be distinguished. 

Figure [2] displays the electrostatic contribution to the 
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excess interfacial tension 7^^ (see Eq. ([T6)) ) of oil droplets 
in water (0/W, ascending curves) and water droplets in 
oil (W/0, descending curves) for an ionic strength in wa- 
ter /w = 1 mM, where U ^ ^^r 0/W and g'f 
for W/0, as a function of the droplet radius a. The 
analytical expression Eq. (I16|) within linearized Poisson- 
Boltzmann theory (thin solid curves) is compared with 
numerical results of the non-linear Poisson-Boltzmann 
theory (thick dotted curves). The slight quantitative dif- 
ferences are due to the linearization approximation and 
they are already present in the planar system (a^^ = 0). 
The quantities n, p, and (pD as well as the crossover val- 
ues i/i and 2/2 correponding to the curves in Fig. [2] are 
displayed in Tab. HI According to Sec. Illli regime III 
is expected to be absent for the W/0 systems because 
Us > > V2 ^ whereas regimes II and III are absent for 
the 0/W systems because > > j/2 (see also the 
discussion of Fig. [S] at the end of this section). 

Due to Eq. pB]) the relative change of the electrostatic 
excess interfacial tension 7|f (a) and the droplet charge 
per droplet surface area a{a) with respect to their pla- 
nar values 7|^(oo) and <t{oo), respectively, are equal, and 
they exhibit the low-curvature asymptotic behavior (see 
Eqs. da and (HIl)) 



G{K\,a,n,p) ,Kba>?/f, (21) 



where G was defined in Eq. ^5)) . Upon rewriting Eq. ([T 
one recognizes the asymptotic behavior 



G{y > y2,n,p) 

I "(I-P') ^-I 
p{l + np) 



n{p + nf _2 
p{l -f npY' 

l + sign(l-n) — 

y 



n{\ — p^) _^ n{p + n)"^ 
p{l + np) ^ p{l + np)'- 



■y 



(22) 



which equals the expansion in y~^ up to second order. 
Hence, the low-curvature expansion up to second order 
in obtained by combining Eqs. (PT|) and is ex- 
pected to be acurate if K\,a ^ j/^ , . Traditionally, em- 
pirically motivated expansions in have been used to 
represent the curvature dependence of the interfacial ten- 
sion without knowing their applicability a priori. How- 
ever, it has been argued by Konig, Roth, and Mecke on 
the basis of a morphometrical approach that the devi- 
ation of intensive thermodynamic quantities from their 
planar values are linear combinations of the mean and 
the Gaussian curvature provided the geometrical length 
scales are much larger than any correlation length [23j . 
i.e., a ^ ^d^'^^b^ with the inverse Debye length in the 
droplet Kd :— K\,p/n, or equivalently n^a ^ yi ,y^ ■ This 
condition is only sufficient but not necessary for the va- 
lidity of the above low-curvature expansion because it al- 
ready implies /tbO ^ 2/2 '^^^ 2/2^ — '^^^(yi^ys)- For 
n,p^ 1, the low-curvature expansion is valid if Kdfl 3> 1, 
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FIG. 3: Electrostatic contribution to the excess interfacial 
tension in mixtures of oil and water (ew = 80) as a func- 
tion of the dielectric constant Eo of the oil for droplet radii 
a G {50 nm, 100 nm, 250 nm, 500 nm, 1000 nm, 00} of an oil 
droplet in water (O/W) and a water droplet in oil (W/0) 
with ion radii a+ = 0.36 nm and a„ = 0.30 nm as well as the 
ionic strength in water Jw = 1 mM. The thin solid curves 
are calculated by means of the analytical expressions within 
the linear theory of Sec. IIIII whereas the thick dashed curves 
are obtained by numerically solving the non-linear Poisson- 
Boltzmann equation ([SJ. There is quantitative agreement for 
all a and Eo considered here. 



independent of the bulk Debye length k.^ ^ , because in this 
case y2 ~ yi ■ This is the case, e.g., for the W/0 systems 
considered in Tab. [H 

From Eq. p8|) one straightforwardly recognizes the 
symmetry G{py/n, 1/n, 1/p) — G{—y,n,p) which means 
that swapping droplet and bulk liquid, i.e., p 1-^ 1/p, 
n '—^ 1/n, Kh I— > Kd, while keeping the droplet radius a 
fixed has numerically the same effect on function G as 
inverting the sign of the droplet radius. Due to this sym- 
metry one concludes for the coefficients of an expansion 
in inverse powers of a as in Eq. (|22p for y — Ki^a that 
upon swapping droplet and bulk liquid the odd-order co- 
efficients merely invert their sign, whereas the even-order 
coefficients do not change. This phenomenon can be ob- 
served in Fig. [21 where the slope close to the planar limit 
(a^^ = 0), which is proportional to the excess Tolman 
length due to the presence of ions, simply changes its 
sign upon swapping oil and water (0/W ^ W/0). 

Figure [3] exhibits the electrostatic contribution to the 
excess interfacial tension 7™ as a function of the di- 
electric constant Eq of the oil for the ionic strength 
in water = 1 mM and for various droplet radii 
a G {50 nm, 100 nm, 250 nm, 500 nm, 1000 nm, 00}. As in 
Fig. [21 the thin solid curves correspond to the analytic 
linear theory of Sec. IIIII whereas the thick dotted curves 
are the numerical results of the non-linear scheme. Quan- 
titative agreement is observed, even in the low-eo range 
where the Donnan potential 4>d is not small and the lin- 
earization approximation is not a priori justified. From 
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FIG. 4; Total excess interfacial tension in mixtures of 
oil and water (sw ~ 80) as a function of the di- 
electric constant So of the oil for droplet radii a G 
{50 nm, 100 nm, 250 nm, 500 nm, 1000 nm, oo} of an oil droplet 
in water (O/W) and a water droplet in oil (W/0) with ion 
radii a+ = 0.36 nm and a_ = 0.30 nm, interfacial width pa- 
rameter |s| = 0.33 nm, as well as the ionic strength in water 
/w = 1 mM. 
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FIG. 5; Electrostatic excess interfacial tension in mixtures 
of oil (eo = 5) and water (ew = 80) as a function of the 
ionic strength in oil (/□) or water (Jw) for droplet radii a G 
{0.1 /im, 1 ^m} of an oil droplet in water (O/W) and a water 
droplet in oil (W/O) with ion radii a+ = 0.36 nm and a_ = 
0.30 nm. The O/W system exhibits only the regimes I and IV 
(see main text and Fig. [l}, whereas for the W/O system the 
regimes I, 11, and IV are present. Upon changing the droplet 
size a the crossover ionic strengths shift by a factor oT^ . 



by 



the linearized theory of Sec. lIIII one can derive the asymp- 
totic behavior = 0{—{eo ~ £w)^) for Eq £w as well 
as = ©(— exp(— const/eo)) for an O/W system and 
71^'= = ©(-I/Eo) for a W/O system as £□ ^ 0. This 
behavior is apparent in Fig. [31 too. 



The total excess interfacial tension 7'''^ comprises not 
only the electrostatic part 7!^ but also the contribution 
7f(f due to the interfacial effects (see Eq. (fTS])). It is 
readily seen that — ±.0{s{e^ — £0)) for £0 — > e„ and 
7i'(f = ±C(s) for £0^0 where the upper {+) and the 
lower (— ) sign correspond to an O/W and a W/O sys- 
tem, respectively. Hence, if s ^ 0, the interfacial effects 
will dominate over the electrostatic effects in the limits 
£0^0 for O/W systems and £0 — > £w for arbitrary 
systems. Figure |4] displays the total excess interfacial 
tension corresponding to the parameters used in Fig. [3] 
and an interfacial width parameter s with \s\ = 0.33 nm 
on the water side of the interface, i.e., s > for O/W 
and s < for W/O. 



According to the results of Sec. IIIII the electrostatic 
excess interfacial tension 7™ as a function of the bulk 
ionic strength /b := g"^^ can be asymptotically described 




(I)/b»/b"i,/, 



b3 



(n)/b"i,/b'2«4«A 



b3 



StT 1 - ?l2 



-alu 



(111)4^, «/b«/, 



b2 



(IV)/b«/b"i 



the 

^2 



crossover bulk ionic 



(23) 



strengths ^y^^ 



with 

^^^y'' ' e (1,2,31, where the y^,k E {1,2,3} are 

defined in Sec. IIIII For a planar system (a = 00) the 
crossovers are at zero ionic strength, hence only the high- 
ionic strength regime I in Fig.[T](/b » ^hi'^hs) present, 
which coincides exactly with the electrostatic contribu- 
tion to the excess interfacial tension in Ref. [S]. 

For an oil dielectric constant £□ = 5 and a droplet 
radius a — 1 /im the crossover bulk ionic strengths are 
^bi « 7lrnM,I^, 
system, where /b 



82nM,4'3 « 93 nM for an O/W 



and I^^ 



/w is the ionic strength in water, 
7.6fM,/^'2 ~ 6.6fM,/^'3 « 5.8 nM for a W/O 



system, where /b = -^o is the ionic strength in oil (see 
Tab.|T|. Here, ionic strengths in oil, /q, and in water, 7^, 
are related to each other by lo/Iw ~ 8.1 • 10~^. Figure[S] 
displays 7™ as a function of the ionic strength in the 
physical range /w S [10"'' M, 10 M] for the droplet radii 
a = 1 /xm and a = 0.1 fj-m. The crossover ionic strengths 
of the latter droplet size are larger by a factor 100 as 



compared to the former because J, 



bfc 



By inspec- 
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tion of the values of the crossover bulk ionic strengths 
one expects only the regimes I and IV of Fig. [T] to be 
present for the 0/W system whereas the regimes I, II, 
and IV are expected for the W/0 system. The occur- 
rence of the regimes I and IV for the 0/W system and 
I, II, and IV for the W/0 system can be inferred from 
Fig. [5] in conjunction with Eq. 



V. CONCLUSIONS AND SUMMARY 

It turned out in the previous section that the ana- 
lytical theory of Sec. IIIII based on a linearized Poisson- 
Boltzmann theory is in good (at least) qualitative agree- 
ment with the results from the full non-linear theory. It 
can therefore be expected that the general conclusions 
drawn from that linear theory apply to more elaborate 
models [13, El, El, too. 

According to Eqs (fT3|) . (|T5l) . and ([T6|) the excess hquid- 
liquid interfacial tension is 7°^^ — 'yf^ + 7^^ where the 
curvature dependence is essentially only due to the elec- 
trostatic part 7°f and not due to the contribution of the 
short-ranged interfacial effects 'jf^. While 7'''^ can in- 
deed be negative, thereby decreasing the total interfacial 
tension, the largest magnitude |7™| is attained at high 
ionic strengths where 7'''^ ~ j^^, i.e., where 7'''^ is essen- 
tially curvature-independent. One has to conclude that 
the unimodal droplet size distribution of W/0 emulsions 
observed by Leunissen et al. cannot be explained by 
the curvature dependence of the interfacial tension due 
to electrostatic effects alone. However, this conclusion 
does not apply to the experiments by Sacanna et al. Q , 
where highly charged colloids instead of monovalent ions 
are present, as the linearized theory of Sec. IIIII is not a 
priori justified for multivalent ions or highly charged col- 
loids. Instead it is an interesting open question to be 



addressed in future studies as to what extent the qual- 
itative low-valency picture drawn here is valid for the 
presence of high- valency particles. 

In summary, the curvature dependence of the elec- 
trolytic liquid-liquid interfacial tension within a simple 
linear Poisson-Boltzmann model in the spherical geome- 
try has been calculated analytically. This linear theory 
turned out to be at least qualitatively reliable as has 
been checked by numerically solving the corresponding 
non-linear Poisson-Boltzmann model. Novel low ionic 
strength regimes, which are not present for a planar 
liquid-liquid interface, have been identified. Low and 
high curvature asymptotics of the interfacial tension have 
been discussed. In particular, it has been found that in 
systems where the ionic strength and the dielectric con- 
stant in the droplet are much larger than in the bulk the 
range of validity of low-curvature expansions up to second 
order in the inverse radius of curvature is independent of 
the bulk Debye length. 
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APPENDIX: SCALING FUNCTION F 

Upon solving the linearized Poisson-Boltzmann equa- 
tion ([U one obtains analytical solutions for the shifted 
electrostatic potential ip, which, via Eq. pH)) . determines 
the scaling function F introduced in Eq. (fTTj): 



where 



F{x,y,n,p) 



Ti 



Ti{p{l + x) hexp 



2yP{l + x))(p{\ + x) + - 
n /V y 



exp 



2/(1 

n 



X. 



pi 



1 r4 + exp(-2y^)r5 



2^^r6 + exp(-2y^)r7 
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^ cosh ^ 
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n J 
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n^[l + x^ h — 
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2\ 



[l + np+ cosh ^ - [n+p + — sinh ^ 

2- 



1 -\- np 



y 



cosh 



n — p 



p{\ — n^) 



,x <Q 



,0; > 0, 



(A.l) 



sinh ^ 



xy\ 
n J 



^ sinh ^ 



n J 



n(l + x) 

py 



(1 - n'^)x 

py 



coshlpxy) + I 1 h X 

py 



n 



{pyf 



smh.{pxy) 
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r5 := n(l + a;) H cosh(pa;y) - 1 + a; H — -^5- smh(pa;y) 

1 + np H j cosh(pa;t/) + \ n■\-p^ j sinh(pa;i/) 

— 1 ) cosh(pa;y) + [n — p ) sinhOpxy). (A. 2) 

Note that F{x, y, n,p) is continuous at x = 0. 
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